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ABSTRACT 

Empirical probability models for BATSE GRB location errors are developed via a Bayesian 
analysis of the separations between BATSE GRB locations and locations obtained with the 
Interplanetary Network (IPN). Models are compared and their parameters estimated using 392 
GRBs with single IPN annuli and 19 GRBs with intersecting IPN annuli. Most of the analysis is 
for the 4Br BATSE catalog; earlier catalogs are also analyzed. The simplest model that provides 
a good representation of the error distribution has 78% of the probability in a 'core' term with 
a systematic error of 1.85 degrees and the remainder in an extended tail with a systematic 
error of 5.1 degrees, implying a 68% confidence radius for bursts with negligible statistical 
uncertainties of 2.2 degrees. There is evidence for a more complicated model in which the error 
distribution depends on the BATSE datatype that was used to obtain the location. Bright bursts 
are typically located using the CONT datatype, and according to the more complicated model, 
the 68% confidence radius for CONT-located bursts with negligible statistical uncertainties is 
2.0 degrees. 

Subject headings: gamma rays: bursts, observations — methods: statistical 



1. INTRODUCTION 

An improved model for the distribution of errors 
for gamma-ray burst (GRB) locations obtained with 
the Burst and Transient Source Experiment (BATSE) 
is presented. The error model should aid in using 
the BATSE locations in projects such as searching 
for counterparts and searching for evidence of repe- 
tition and clustering, including gravitational lensing. 
Most of the analysis applies to the 4Br catalog (Pa- 
ciesas et al. 1998), which has some revised locations 
compared to earlier catalogs, including the initial 4B 
catalog (Paciesas et al. 1997), which was released 
electronically and on CDROM. Summary results are 
presented for earlier catalogs. 

The BATSE instrument consists of eight modules 
located on the corners of the Compton Gamma-Ray 



Observatory (CGRO). Each module contains a Large 
Area Detector (LAD), consisting of a 1.27 cm thick 
by 50.8 cm diameter Nal crystal viewed by three 
photomultiplier tubes. Each LAD has an approxi- 
mately cosine response. GRB locations are deter- 
mined with the program LOCBURST by modeling 
the responses of either four or six detectors as a func- 
tion of assumed source location, intensity and spec- 
trum (Pendleton et al. 1998). The detector response 
model is based upon Monte Carlo photon propaga- 
tion, laboratory measurements and space observa- 
tions. It includes scattering from nearby spacecraft 
structures and from the earth's atmosphere (Pendle- 
ton et al. 1995, 1998). 

The error models presented in this paper are 
obtained empirically by statistical comparison of 



'Department of Physics, University of Alabama in Huntsville, Huntsville, AL 35899 

2 Center for Space Plasma, Aeronomic and Astrophysics Research, University of Alabama in Huntsville, Huntsville, AL 35899 
3 Space Sciences Laboratory, University of California, Berkeley, CA 94720-7450 
4 NRC & NASA/Marshall Space Flight Center, Huntsville, AL 35812 
5 NASA/Marshall Space Flight Center, Huntsville, AL 35812 



2 



BATSE GRB Location Errors 



BATSE locations with locations obtained by other 
techniques. The ideal reference dataset would have 
many locations to enable a good statistical compar- 
ison and would consist of locations with errors sig- 
nificantly smaller than the BATSE errors so that all 
discrepancies could be attributed to BATSE. If such 
a dataset existed, there would probably be little inter- 
est in the BATSE locations. The reference set we use 
is the locations of the Interplanetary Network (IPN) 
Supplements to the BATSE 4Br catalog (Hurley et 
al. 1998a,b; Laros et al. 1997,1998; Hurley 1998). 
Although this dataset is large, most of the IPN loca- 
tions are highly accurate only in one dimension. 

2. ANALYSIS APPROACH 

2.1. The Problem 

The goal is an improved model of the error distri- 
bution of BATSE GRB locations. The error model 
is a probability density function p defined on the 
sphere. A probability density function can have any 
nonnegative value; integrating p with respect to solid 
angle element dVL over any solid angle region VL 
yields the probability P (between and 1) that the 
true location is in the region Q. 

A portion of the errors originates from the Poisson 
fluctuations of the detected counts. These fluctua- 
tions are propagated into an estimated statistical un- 
certainty a stat by the program LOCBURST. The loca- 
tions are found by LOCBURST by x 2 -minimization, 
and cr stat is calculated as the radius of the spherical 
small circle with the same area as the 68% confi- 
dence ellipse obtained from the Hessian matrix of 
derivatives of x 2 (Pendleton et al. 1998). In this 
paper we assume that the total error distribution is 
azimuthally symmetric. The value of cr sta t for each 
GRB is listed in the 4Br catalog (Paciesas et al. 
1998), or, for bursts with unchanged locations, in the 
3B catalog (Meegan et al. 1996). We attribute the 
difference between the total error cr tot and the statis- 
tical uncertainty cr stat to a systematic error cr sys . Pos- 
sible sources of systematic errors include: inaccura- 
cies in the assignment of energies to channel bound- 
aries, inaccuracies in the deadtime correction of high 
count rates, the approximations of the statistical and 
total error boxes as circles, inaccuracies in the detec- 
tor response model as a function of photon energy 
or direction, including inaccuracies in the model of 
scattering from the spacecraft and the earth's atmo- 
sphere, deviations of the actual GRB spectra from the 
power law assumed in most cases, and inaccuracies 
in background subtraction. 

The first published BATSE catalog, the IB catalog 



(Fishman et al. 1994), had 260 GRBs. The error 
estimate was based on eleven events with small error 
boxes from IPN data. A root-mean-square difference 
between the actual BATSE-IPN separations and the 
BATSE statistical uncertainties estimated a system- 
atic error of a sys = 4° (Fishman et al. 1994). The 
3B catalog (Meegan et al. 1996) introduced an im- 
proved location algorithm. The analysis of the 3B 
location errors expanded the reference dataset to in- 
clude 38 locations obtained with CGRO/COMPTEL 
or WATCH in addition to 12 IPN locations. The 
systematic error cr sys was estimated to be 1.6° from 
the root-mean-square difference between the actual 
separations to these 50 reference locations and the 
BATSE statistical uncertainties. The total errors of 
the COMPTEL or WATCH locations are not much 
smaller than the BATSE errors; in this situation a 
very good understanding of the errors of the refer- 
ence locations is required, otherwise the portion of 
the location discrepancies due to BATSE may be 
miscalculated. An overestimation of the COMP- 
TEL or WATCH location errors may have caused 
the BATSE systematic error to be somewhat under- 
estimated. In addition, we show that the BATSE er- 
ror distribution is not well characterized by a single- 
valued systematic error so that various techniques of 
estimating a single systematic error yield different 
values (see, e.g., Pendleton et al. 1998). 

We term the simple error model specified in the 3B 
paper (Meegan et al. 1996), that the total error cr tot 
is formed by adding the statistical uncertainty cr stat 
and a systematic error a sys of 1.6° in quadrature, the 
"minimal" model. 

A rough rule of thumb learned in this analysis is 
that with /V reference locations, models with up to 
~ N/50 parameters may be fit. Exceeding this em- 
pirical limit typically results in models difficult or 
impossible to fit, and, if a fit is obtained, nonsense 
parameter values. There are only about 20 'point' lo- 
cations available from IPN data, but by using a more 
sophisticated analysis method (Graziani and Lamb 
1996), the several hundred single annuli in the IPN 
Supplements may be used to constrain the BATSE 
location error distribution. Formerly representations 
of the systematic error were limited to using a sin- 
gle parameter; considerably more complicated error 
models can now be tested. 

2.2. Reference Data 

The IPN location technique compares the arrival 
times of the GRB wavefront at widely separated in- 
struments, thereby determining the angles the plane 



Briggs et al. 



3 



wave makes with the vectors connecting each pair of 
instruments (Hurley et al. 1998a). With only one 
pair of spacecraft, one angle is determined, and the 
error region is an annulus. This is the typical case 
during the time period of the 4B catalog. Ideally, 
three or more spacecraft have separations of inter- 
planetary scale, multiple annuli are found, and their 
intersections give a small error box. The IPN Supple- 
ments to the 4Br catalog contains 458 annuli for 412 
GRBs. Of these 412 GRBs, 366 were observed by 
only one deep- space instrument, Ulysses (Hurley et 
al. 1998a,b). The remaining 46 GRBs were observed 
by Ulysses and another deep-space instrument such 
as Pioneer Venus Orbiter or Mars Observer (Laros 
et al. 1997,1998). 

Measurements of the actual separations 7 are avail- 
able for the few cases with small error boxes from 
intersecting IPN annuli. For the cases for which 
only single annuli are available, the BATSE-IPN sep- 
arations are characterized by the distances of clos- 
est approach p of the annuli to the BATSE locations 
(Fig. 1). 

Our analysis approach assumes that the reference 
locations are effectively points or lines compared to 
the BATSE locations, an assumption which is ful- 
filled by eliminating all annuli with 3 a widths ex- 
ceeding 1.6° and eliminating 3 a error boxes with any 
corner-to-corner dimension exceeding 1.6°. Because 
there are few IPN annuli with large widths, decreas- 
ing or increasing the 1 .6° requirement by a factor of 
two leads to only small changes in the number of an- 
nuli included in the comparison sample and negligi- 
ble changes in the results. We convert the annuli to 
circles on the sphere, or their intersections to points, 
by using the circle in the middle of each annulus. 
Each pair of annuli has two intersections — we as- 
sume that the intersection closest to the BATSE loca- 
tion is the true location. Many of the intersecting an- 
nuli create long error boxes which are rejected by the 
requirement that no dimension exceed 1.6°, in which 
case we keep only the narrowest annulus. Such er- 
ror boxes are common during the time period of the 
4Br catalog because frequently several of the space- 
craft of the Third Interplanetary Network have been 
close in space. The vectors connecting each pair of 
spacecraft define the centers of the annuli, and with 
some spacecraft in similar directions, these vectors 
are typically close to parallel. This leads to 'graz- 
ing' intersections of the annuli and long error boxes. 
(An additional consequence of grazing intersections 
is that a small error in one annulus will move the 
intersection a comparatively large distance). After 



removing the wide annuli and large error boxes, 430 
annuli for 411 GRBs remain. For GRBs with inter- 
secting annuli meeting our criteria and therefore with 
a measurement of 7, we do not use the two measure- 
ments of p that could be obtained from the individual 
annuli, thereby ensuring the independence of all the 
7 and p measurements. These criteria yield a ref- 
erence dataset with 19 measurements of 7 and 392 
measurements of p. 

The 411 GRBs with IPN locations meeting our cri- 
teria represent 25% of the 1637 GRBs of the 4Br 
catalog. The relative sensitivities of BATSE and the 
deep- space instruments of the Third Interplanetary 
Network are largely determined by their effective ar- 
eas and by their backgrounds. Each BATSE LAD 
has an effective area of 2025 cm 2 , compared to, for 
example, a projected area of 20 cm 2 for Ulysses. 
The differences in sensitivities are not as great as 
the comparisons of areas would indicate because the 
deep- space instruments have lower and more stable 
backgrounds. Fig. 2 compares the properties of the 
bursts of the 4Br catalog with the bursts of the IPN 
Supplements to the 4Br catalog. The IPN Supple- 
ments sample almost the entire range of peak fluxes 
on the 64 ms timescale of the 4Br catalog, albeit 
with emphasis on the bright bursts. In terms of flu- 
ence, the IPN Supplements have a strong bias to- 
wards bright bursts, with only rare events below 10" 6 
erg cm" 2 having IPN annuli. There is a close corre- 
lation between burst fluence and BATSE cr stat , so a 
similar effect is seen in the distributions of cr stat . 

2.3. Bayesian Model Comparison 

If many 'point' reference locations were avail- 
able, a histogram of the observed separations 7 be- 
tween the BATSE locations and the reference lo- 
cations would provide a good representation of the 
probability distribution of BATSE errors. Because 
most of the reference locations from the IPN Sup- 
plements are single annuli, insufficient information 
exists for such a simple approach. Instead the error 
function 77(7) is determined by comparing models for 
p(j) and optimizing their parameters using the entire 
dataset of 7 and p measurements. 

A Bayesian analysis very similar to that of 
Graziani and Lamb (1996) is used to identify the 
best error model. The Bayesian approach identifies 
the likelihood function as the correct 'merit' function 
(Sivia 1996). The likelihood function L is simply the 
product of the probability density functions p eval- 
uated at each observed value of 7 for bursts i with 
intersecting annuli and at each observed value of p 



4 



BATSE GRB Location Errors 



T 



Fig. 1 . — If two thin IPN annuli are available (solid and dot- 
ted arcs), then we assume that their intersection gives the true 
location T. In this case the angle 7 between the true location 
T and the BATSE location B is known. If only one IPN annu- 
lus is available (solid arc), the position of the true location on 
the annulus is unknown, and the separation between the annulus 
and BATSE location B is characterized by the angle of closest 
approach p. 



for bursts j with single annuli: 

L = I[Pi(li)I[Pp(Pj)- 



(1) 



Combining these considerations, the Odds ratio O b /a 
specifies the factor by which one should prefer model 
B over model A: 



The probability density function p p used for the sin- 
gle annuli cases is analytically derived from p 7 (see 
Appendix). The likelihood function identified by the 
Bayesian method specifies how measurements of dif- 
fering nature (intersecting vs. single annuli) are to 
be combined in the analysis. While eq. 1 shows the 
likelihood L as a function of the observations, L is 
also, through the functions p 7 and p p , a function of 
the error model and its parameters. For a particular 
error model, the best-fit parameters are obtained by 
maximizing the likelihood. 

An additional consideration arises in comparing 
models: a model with more parameters than another 
must be penalized because it has greater ability to 
match the data whether or not it is correct. The 
Bayesian approach specifies how this must be done: 
each model is penalized by an Occam's factor F. As- 
suming that prior to examining the data one expects 
parameter X k to be in the range A™" to Xf 3 * and that 
the likelihood function is approximately Gaussian, 
these factors are, per parameter X k (Sivia 1996), 



Ob/a = 



P(B) /> ior (B) x L(B) x F(B) 



F k = 



\ max \ min ' 



(2) 



where a\ k is the uncertainty on A^ obtained from the 
fit. The overall Occam's factor F is the product of the 
F k and penalizes a model for obtaining a better fit by 
using adjustable parameters, with a greater penalty 
(smaller F) when the parameters have larger possi- 
ble prior ranges. 

The final consideration is that one might have prior 
preferences for particular models from earlier analy- 
ses or judgments . These preferences can be repre- 
sented as prior probabilities for each model P piior (M). 



P(A) P prior (A)xL(A)xF(A)' 



(3) 



To avoid suspicions that our beliefs have created the 
results, we equate all prior probabilities so that the 
first factor in eq. 3 cancels. Many of the models 
have common parameters so that some of the Oc- 
cam's factors F k also cancel. For most of the model 
comparisons the likelihood ratios are so large that no 
reasonable priors nor Occam's factors F could alter 
which model is identified as best. 

Further discussion of Bayesian methodology, 
model fitting and comparison are given by Sivia 
(1996) and Loredo (1990). 

Instead of listing the odds ratios b /a for every 
combination of models A and B, we list what we call 
odds factors M for each model M. The odds factors 
M are proportional to P(M) by a factor which is the 
same for each model, therefore the odds ratios Ob/a 
can be found by forming the ratio of the factors Ob 
and A , or by differencing their logarithms (Tables 
2 and 5). The likelihoods and parameter uncertain- 
ties are also listed (Tables 2 and 3) so that others can 
calculate odds ratios based upon their priors. 

While the quantitative analysis is done with un- 
binned data using likelihood, we prefer to depict the 
data and models with histograms (Fig. 3-5, 7 & 9). 
In order to present the results in a consistent manner, 
the data and models are binned using the total errors 
of the minimal model, cr° t = (a s 2 tat + 1.6° 2 ) 1/2 . 



2.4. Probability Models 
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FlG. 2. — Comparison of bursts for which BATSE triggered (solid histograms) with the subset in the IPN Supplements (dashed 
histograms), a) The distributions of peak flux as observed by BATSE between 50 and 300 keV on the 64 ms timescale. There are 
1292 bursts in the 4Br catalog with flux/fluence values (solid histogram), of which 324 are also in the IPN Supplements (dashed 
histogram), b) The distributions of fluence as observed by BATSE between 50 and 300 keV of the same groups of 1292 and 324 
bursts, c) The distributions of BATSE er stat for all 1637 GRBs of the 4Br catalog (solid histogram) and for the subset of 41 1 events 
with IPN data meeting our criteria (dashed histogram). 
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We define the total error as the quadratic sum of 
the statistical uncertainty cr stat (as listed in the 4Br 
catalog) and a systematic error: 



„ (JL ,J2 \l/2 
<7 tot = (^ st at + ^ S ys) 



(4) 



Our error models are based on the Fisher proba- 
bility density function p F , which has been called the 
Gaussian distribution on the sphere: 



7? F (7) dVL = 



K 



27r(e K -e K ) 



3 KCOS7 



(5) 



where 7 is the angle between the measured and true 
location (Fisher et al. 1987). The first convenience 
of this distribution is that it is normalized on the 
sphere and that the 'volume' element is clearly solid 
angle dVl. The distribution also has convenient ana- 
lytic properties; in the small angle approximation it 
reduces to the Gaussian distribution. An algorithm to 
simulate random locations from a Fisher distribution 
is given by Fisher et al. (1987). 

The probability P of the true location lying in a 
region of solid angle Q is found by integrating the 
probability density function p: 



P = 



K 



27r(e K -e K ) Jn 



cos 7 



(6) 



Specializing to the probability of the true location be- 
ing in the ring 71 < 7 < 72, eq. 6 becomes 



P = 



/ d7sin 7 e KCOS7 . 

K — e K J-yi 



(7) 



The parameter k is termed the concentration pa- 
rameter. BATSE GRB location errors have tradition- 
ally been specified in terms of a, which has been de- 
fined by the BATSE team as the radius of the circle 
with the same area as the 68% confidence ellipse. 
Assuming the error box to be circular, and setting 
P = 0.68, 71 = and 72 = cr tot in eq. 7, one obtains, for 
cr tot in radians and for er tot small (<; 20°): 



1 



K = 



(0.66a tot ) 2 



(8) 



Models more complicated than a Fisher distribu- 
tion are tested. The first method of building more 
complicated models is to make cr sys a function of ei- 
ther intrinsic burst properties or of instrumental fac- 
tors which might influence the quality of locations. 
The second way in which more complicated models 



are formed is to sum two Fisher distributions with 
different values of cr s 



'sys- 



P = /iPf + (1-/i)Pf- 



(9) 



In this case, for each burst, the values of the two total 
location errors <7 t0ti , are calculated from the statistical 
uncertainties cr stat listed in the catalog and the model 
values cr syS)i (for i =1,2): 



Ctot,/ - (°"stat + ^sys,/) 



A/2 



(10) 



For the single annulus cases, the probability den- 
sity function p(p) is analytically derived from the 
probability density function ^(7) (see Appendix). 
Since p < 7 (see Fig. 1), P(p < X) > P{^ < X). Ta- 
ble 1 gives some cumulative probabilities P for both 
Pvin) and the corresponding p(p). 

3. RESULTS 

3.1. 4Br Catalog 

The systematic error of 1.6° used in the minimal 
model was determined using 50 comparison loca- 
tions, some of which had errors not much smaller 
than the BATSE errors (Meegan et al. 1996). Fig. 3 
compares the 4Br data to the minimal model — with 
the large and accurate IPN reference dataset, the min- 
imal model is clearly seen to be a poor fit to the 
data. Model 1 generalizes the minimal model by op- 
timizing the value of the systematic error, obtaining 
2.8°. While a major improvement over the minimal 
model, this model fits the data poorly (Fig. 4) and has 
a much smaller odds factor than some other models 
(Table 2). 

The value of 1.6° published in the 3B catalog 
(Meegan et al. 1996) was determined by a root- 
mean-square comparison of the observed separations 
with the combination of the BATSE statistical uncer- 
tainties and the known errors of the reference loca- 
tions. The 3B catalog warned that some locations 
might be substantially worse than average and that 
the error distribution might have a non-Gaussian tail. 

Model 2 implements the idea of a Gaussian-like 
distribution with an extended tail by summing two 
Fisher distributions having differing systematic er- 
rors (eqs. 9 and 10). This model provides an excel- 
lent fit to the data (Fig. 5 and Table 2). The evidence 
favors Model 2 over the minimal model by an odds 
ratio of 10 56 and over Model 1 by 5 x 10 7 (Table 2). 
The best- fit parameters are given in Table 3. 

Model 2 is best thought of as a 'core-plus-tail' rep- 
resentation of the error distribution rather than as a 
two-component model. The model does not specify 
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Table 1 

Cumulative Values for the Distributions of 7 and p assuming the Fisher Distribution 



Angle X P F ( 7 < X) P(p < X) 



\a 0.6826 0.8703 

2a 0.9898 0.9976 

3a >0.9999 >0.9999 




FlG. 3. — Histograms of the data (solid) and Model (dashed) for the separations 7 between BATSE locations and intersecting 
IPN annuli and the closest separations p between BATSE locations and single IPN annuli. The data and model are binned in terms 
°f °tot = ( <J s 2 tat+ 1 -6° 2 ) 1//2 - The data and Model are clearly in poor agreement. 
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Fig. 4. — Histograms of the data (solid) and Model 1 (dashed). Model 1 predicts too few events at both small {p/a^ <J 0.25) and 
large (p/<r t " t ^ 3) separations. 



o 1 



y / °tc 




Fig. 5. — Histograms of the data (solid) and Model 2 (dashed). The agreement between the data and Model 2 is excellent. 
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Table 2 

Model Comparison Results for 4Br Catalog 



Model 


Description 


Number of 
Parameters 


likelihood 


togio 
Odds Factor 





Minimal: Systematic 
Error Assumed to be 1.6° 





161.5 


161.5 


1 


Single Systematic Error 
applied to all locations 


1 


211.5 


210.2 


2 


Core-plus-Tail Systematic 
Errors 


3 


220.8 


217.9 


4 


Systematic Error a power 
law function of statistical 
uncertainty: 

c%s =Ao stat /i°r 


2 


212.6 


210.4 


10 


Core-plus-tail Systematic 
Errors depending on 
datatype 


6 


223.7 


219.4 


12 


Core-plus-tail for CONT, 
one systematic error for 
other datatypes 


4 


222.9 


219.0 


33 


Core-plus-tail, only core 
systematic error depends 
on datatype 


4 


223.6 


219.9 
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Table 3 

Best-Fit Error Model Parameter Values for the 4Br Catalog 



Model 


Symbol 


Parameters 
Description 


Valued 


1 


<^sys 


Systematic Error 


2.80±0.12 


2 


^sys 
r 
fl 

a 2 

sys 


Core systematic error 
rraction in Core term 
Tail systematic error 


1.85° ±0.16° 

A TOiA AO 

0.78 0.08 

c 1 o +0.8 
D - L -0.6 


4 


A 

a 


Amplitude (error for cr stat = 1°) 
Power law index 


2.93±0.13 
0.14 ±0.06 


10 


°CONT 

/cONT 

2 

a CONT 
1 

° other 
f l 

J other 
2 

°other 


Core systematic error for CONT 
Fraction in Core Term for CONT 
Tail systematic error for CONT 

Core systematic error for non-CONT 
Fraction in Core Term for non-CONT 

Tail systematic error for non-CONT 


1.68° ±0.16° 
0.82881 

C OO +1.0 

J - J -0.8 

2.6° ±0.5° 

A 7 o+0.19 

c 00+2.1 
D - Z -1.2 


12 


°CONT 

/cONT 

2 

°"C0NT 
°other 


Core systematic error for CONT 
Fraction in Core Term for CONT 
Tail systematic error for CONT 
Systematic error for non-CONT 


1.67° ±0.15° 
0.82881 

C OO +1.0 

J - J -0.8 

3.44° ±0.26° 


33 


_1 

CONT 

^other 

/all 
7. 

^all 


Core systematic error for CONT 
Core systematic error for other 
Fraction in Core, all datatypes 
Tail systematic error for all datatypes 


1.67° ±0.15° 
2.74° ±0.34° 

0.82885 

5.4° ±0.8° 



tThe uncertainties on the parameter values are for single parameters of interest (i.e., they were obtained from 
the change in likelihood equivalent to the usual x 2 + 1 prescription). In some cases the uncertainties between 
the parameters are highly correlated. However, further specification of the errors on the errors of the locations 
would be excessive. 
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that some bursts belong to a core component (small 
cr sys ) and others to a tail component (large cr sys ); in- 
stead all bursts have a tail to their location error dis- 
tribution. Model 2 uses no burst property other than 
cr stat to determine the location error distribution of 
each burst. 

To depict the distribution of separations 7 between 
the BATSE location and the true location given by 
Model 2, we perform part of the solid angle integra- 
tion over the sphere of the probability density func- 
tion p 2 of Model 2: 

1= / dj difj sin^ p 2 in) (11) 
Jo Jo 




The function dP 2 /d r y is depicted in Fig. 6 for a burst 
with a negligible value for cr stat . The 'core' and 
'tail' terms are shown separately along with the to- 
tal model. While 22% of the probability is in the tail 
term, much of the area of the tail term is near 7 = 
so that only 7% of the locations are past 7 = 5°. 

We have endeavored to find correlations that would 
allow us to assign a sys values based upon burst or lo- 
cation properties, either by fitting models with a sys 
a function of the property, or by dividing the loca- 
tions into categories by the property and comparing 
the values of model parameters. Among the proper- 
ties for which no significant dependence was found 
are: the fluence of the burst, the value of <7 stat , the x 2 - 
value obtained by LOCBURST, the angle to the cen- 
ter of the earth (which affects the importance of at- 
mospheric scattering), the spectral index of the burst, 
and whether the event was a Ulysses trigger or was 
obtained from Ulysses continuous data. In a test for 
a dependence of location quality on spacecraft scat- 
tering, we made three extensions each of Models 1 
and 2, adding a parameter dependence on the space- 
craft hemisphere (+/-X, +/— Y, or +/-Z) in which 
the burst was located. The coordinate system is de- 
fined so that COMPTEL and EGRET view in the +Z 
direction and OSSE scans in the XZ plane (Gehrels, 
Chipman & Kniffen 1993). The motivation is that 
the spacecraft is not symmetrical in Y and Z, and in- 
adequacies in the modeling of the scattering from the 
spacecraft might cause location quality to depend on 
position in spacecraft coordinates. There were mod- 
est odds ratios improvements for two of the exten- 
sions of Model 1 (not approaching the odds ratio of 
Model 2), but none for the extensions of Model 2. 
We conclude that the bilateral asymmetries of the 



spacecraft do not create bilateral asymmetries in lo- 
cation quality. 

The only property that we have identified as sig- 
nificant is the type of the data used to obtain the 
location. In addition to four continuously transmit- 
ted 'background' datatypes, the BATSE instrument 
packages data from triggered events into eight burst 
datatypes. Any of the background or burst datatypes 
which include data separately from each of the LADs 
can be used to derive a location. The default datatype 
for locating bursts is CONT, which provides contin- 
uous coverage regardless of trigger status in 16 en- 
ergy channels with 2 s time resolution. The CONT 
datatype was used for 58% of the locations in the 
4Br catalog. If the event is shorter than 2 s, so that 
using CONT data would unnecessarily add statistical 
fluctuations, or if CONT data are unavailable due to 
a telemetry gap, discriminator data (DISCLA, DIS- 
CLB, PREB or TTE) with four energy channels are 
used. If normal telemetry is available, 1 s resolu- 
tion DISCLA would typically be used. A duplicate 
of DISCLA data, DISCLB, is included in the burst 
data to obtain more reliable telemetry. Together DIS- 
CLA and DISCLB data were used for 23% of the 
locations. The instrument triggers at the end of the 
interval in which a significant rate increase is regis- 
tered. If the burst has essentially ended by the trig- 
ger time, then one of the high-time resolution pre- 
burst datatypes is used, either PREB or TTE. This 
was the case for 12% of the locations. Finally, if no 
regular data are available due to missing telemetry, 
the maximum rates recorded by the instrument dur- 
ing the burst readout are used (MAXBC or MAXC1). 
This was necessary for 8% of the 4B locations. (The 
percentages sum to 101% due to rounding.) 

Models 10, 12 (§3.2) and 33 incorporate cor- 
relations between location accuracy and datatype. 
Model 10 is very similar to Model 2, with the differ- 
ence that the model parameter values depend on the 
datatype. One set of model parameters like those of 
Model 2 are used for locations obtained using CONT 
data, cr^oNT' /cont' an d 0cont> while another set, 
^dther' /other' an d °"other' * s use d for all other locations. 
Model 10 has a likelihood x800 larger than that of 
Model 2 (Table 2). Considering the three additional 
parameters, the odds ratio favors Model 10 by a fac- 
tor of 30. 

Two pairs of parameters of Model 10, flom & 
f^ thei and crc ONT & cither' nave values which are con- 
sistent (see Table 3), i.e., the data do not demonstrate 
that these parameters have different values. This sug- 
gests a model with fewer parameters: Model 33 has 
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Fig. 6.— The function dP 2 /dj (eqs. 11 and 12) is the er- 
ror distribution versus separation 7 according to Model 2. It 
is shown for a gamma-ray burst with a negligible value for <7 stat . 
The function as depicted includes the solid angle factor 27rsin?/>. 
If this factor were not included, the function would be very large 
near 7 = 0, yet the probability of very small 7 is low because of 
the solid angle factor. The integral from 7 = to 180° of the 
function shown is one. The solid curve shows Model 2 (as spec- 
ified in Table 3), while the dashed curves show the two terms, 
'core' and 'tail', separately. 



y (degrees) 



distinct parameter values for the core systematic er- 
rors of CONT-derived locations and locations ob- 
tained with other data types, but only a single param- 
eter fl xx representing the probability of the core term, 
and only a single tail systematic error cr^. The likeli- 
hood ratios of Models 10 and 33 are virtually identi- 
cal; because of the fewer parameters of Model 33 its 
odds ratio is x3 better (Table 2). This improvement 
may mean that the 'tail' is caused by factors indepen- 
dent of datatype, or simply that there is insufficient 
data to discern differences in the tail as a function 
of datatype. The odds ratio improvement is marginal 
and the models make very similar predictions, so an 
explanation is not required. Model 33 is favored over 
Model 2 by an odds ratio of 100. Qualitatively, this 
might be termed as persuasive but not compelling ev- 
idence for Model 33 over Model 2. The data and 
Model 33 are compared in Fig. 7. 

Locations obtained with datatypes other than 
CONT, i.e., using either PREB, TTE, DISCLA, DIS- 
CLB, or MAXBC datatypes, are consistent with be- 
ing of equal quality. However, of the 41 1 bursts with 
reference IPN locations, only 110 were located with 
datatypes other than CONT, so the comparisons are 
probably insensitive. 

The Bayesian approach relies upon model com- 
parisons rather than goodness-of-fit tests — generally 
one does not know whether a particular best-fit value 
of the likelihood is reasonable. As a final test of 
the quality of the models, we (not being doctrinaire 
Bayesians) perform goodness-of-fit tests via simula- 



tions. Assuming a specific model to be true, we sim- 
ulate datasets of BATSE locations using the best-fit 
parameters of the model. We use the actual IPN an- 
nuli, pick locations on single annuli and on the in- 
tersection of each pair to be the 'true' locations, and 
then create simulated BATSE locations using the as- 
sumed error model. The separations between the an- 
nuli and the simulated locations yield simulated val- 
ues of p and 7. The simulated dataset is fit with the 
assumed model, re-optimizing its parameters and ob- 
taining the best-fit value of the likelihood. This pro- 
cess is repeated to generate 1000 simulated datasets 
and likelihood values. If the actual best- fit likelihood 
value is much different from the simulated best-fit 
likelihoods, the model is probably a poor explana- 
tion of the data. The simulation results are shown 
in Table 4 — all of the actual best-fit likelihoods are 
reasonable according to the simulations except for 
the likelihood of the minimal model. Even Model 1 
passes this goodness-of-fit test despite the strong re- 
jection of the model by the Bayesian model compar- 
isons. 

Recommendations. Models 2 and 33 both make 
excellent fits to the data (Tables 2 and 4) and the two 
model histograms are extremely similar (Fig. 5 and 
7). While there is evidence that Model 33 is bet- 
ter, for most applications we recommend Model 2 
because of its simplicity. For critical analysis, 
Model 33 might be selected. The models are speci- 
fied by eqs. 4 to 10. The best-fit parameter values are 
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FlG. 7. — Histograms of the data (solid) and Model 33 (dashed). While quantitative measures show Model 33 to be better than 
Model 2, the model histograms are almost identical. 



Table 4 

goodness-of-flt tests via simulations 







Mean of 


Standard 


Fraction of 




Best Fit 


Simulations 


Deviation 


Simulations with 


Model 


lo gio L fit 


loglO^simu 


of log 10 L simu 


^simu ^ Lfit 





161.5 


293.8 


6.5 


0/1000 


1 


211.5 


214.9 


6.3 


0.29 


2 


220.8 


221.6 


9.0 


0.46 


4 


212.6 


213.4 


6.2 


0.45 


10 


223.7 


222.2 


8.5 


0.58 


12 


222.9 


221.0 


8.2 


0.60 


33 


223.6 


222.8 


8.8 


0.54 
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listed in Table 3. If only confidence radii are needed, 
these may be read from Fig. 8. Because the models 
are empirically determined using 411 reference lo- 
cations, confidence radii much beyond 99% have the 
character of an extrapolation and should be used with 
caution. 

3.2. Comparisons with Previous Analyses 

In a preliminary version of this project (Briggs 
et al. 1998), similar results were obtained. The 
largest difference is that Model 10 was not presented 
in that work because it was inferior to Model 12 and 
Model 33 was not presented because it had not yet 
been tried. Model 10 has two values of <r sys depend- 
ing on datatype (6 parameters total), while Model 12 
has two values of <7 sys for CONT-based locations, but 
only one value for non-CONT locations (4 param- 
eters total). Model 33 has two values for the core 
systematic error, but a single value for the fraction 
and a single-value for the tail systematic error. The 
current analysis (Table 2) shows Model 33 to be bet- 
ter than Model 12 by a small odds ratio of 8. At this 
level additional data or different priors might cause 
the model identified as best to change-this also ap- 
plies to closely related models not presented. This is 
not a problem because these three models are based 
upon the same fundamental ideas (core-plus-tail dis- 
tribution and datatype dependence) and make nearly 
identical predictions for location errors. 

Reasons to prefer Model 10 are that it seems very 
plausible that the error distribution of both CONT 
and non-CONT locations should be best-fit with a 
core-plus -tail model and there is no reason that the 
fractions and tail systematic errors must be identical. 
(This assertion of reasonableness might be termed a 
"prior", but we chose not to quantify our inclination.) 
However, we have recommended Model 33 because 
it has the best odds ratio and is simpler than Mod- 
els 10 and 12; in any case these three models provide 
very similar predictions. 

The parameter values of the models are also 
slightly changed from the preliminary analysis. The 
differences between the analyses are due to improve- 
ments in the IPN locations between the preliminary 
data releases and the published values, improve- 
ments in the BATSE locations between the 4B and 
4Br catalog (Paciesas et al. 1998) (also see §3.3) 
and better rejection of long IPN error boxes. The re- 
moval of a small number of outliers due to mistakes 
had a disproportionate effect on the error analysis. 

An earlier work (Graziani and Lamb 1996), which 
introduced the analysis method used herein, reached 



substantially different conclusions. Their best mod- 
els were intensity dependent, with <r sys a power-law 
function of <7 sta t (Model 4) or of the fluence of the 
burst (Model 8). Table 2 shows Model 2 to be supe- 
rior to Model 4 by an odds ratio of 3 x 10 7 and the 
histograms comparing the observations and Model 4 
(Fig. 9) show the fit to be poor. The fluence-based 
model cannot be analyzed for the full 4Br catalog be- 
cause fluence values are unavailable for some bursts 
due to telemetry gaps that are uncorrelated with burst 
properties. Of the 411 bursts of our comparison sam- 
ple, 324 have fluence measurements. Table 5 lists 
model comparisons using these 324 GRBs. Model 2 
is superior to Model 8 by an odds ratio of 1 x 10 8 . 
The parameter values found for Model 8 are A = 
2.66° ±0.12° and a = -0.084 ±0.036. Like Graziani 
and Lamb (1996), we find the intensity-dependent 
models to be favored over Model 1 (Tables 2 and 5). 

Graziani and Lamb (1996) analyzed the 3B cat- 
alog (Meegan et al. 1996) (the 4Br catalog did 
not yet exist) and restricted the analysis to events 
with single annuli. They neglected all of the events 
with 7 measurements in order to avoid any locations 
that might have been used to optimize the location 
program LOCBURST. Twelve GRBs with 7-values 
were used to identify which proposed algorithmic 
improvements were worth implementing (Pendleton 
et al. 1998). All location algorithm parameter values 
were determined from laboratory measurements, in- 
orbit observations of the Crab Nebula and pulsar and 
from Monte Carlo simulations-none were optimized 
using GRB measurements. We know of no reason to 
exclude all 7 measurements and judge that there is 
little 'circularity' in using the 7 measurements of the 
twelve 'study' bursts. In any case, excluding the data 
for these twelve events has negligible effects on the 
results. 

The reasons that are probably most important in 
explaining the difference between our results and 
those of Graziani and Lamb (1996) are the mod- 
els tested and revisions to the preliminary IPN data 
(Hurley et al. 1998a). Using the original 3B BATSE 
data and the data selections of Graziani and Lamb 
(i.e., excluding bursts with 7 measurements), but us- 
ing the finalized IPN data, Model 2 is favored over 
Model 4 by an odds ratio of 400. Adding the 7 mea- 
surements for all except the twelve 'study' bursts, 
the odds ratio becomes 4700. Unlike Graziani and 
Lamb, we conclude that there is no evidence for a 
direct dependence of BATSE systematic location er- 
rors on burst intensity. 

Because the IPN Supplement emphasizes rela- 
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FlG. 8. — Confidence Radii for GRB location as a function of cr stat listed in the 4Br catalog. The curves are based upon Models 2 
and 33, as specified in Table 3. The solid curves show the confidence radii for Model 2, the dotted curves the confidence radii 
according to Model 33 for locations obtained using the CONT datatype, and the dashed curves the radii according to Model 33 for 
locations obtained using datatypes other than CONT. The bottom trio of curves shows the radii in which the location will fall 68% 
of the time according to the models, the middle trio the 90% confidence radii, and the top trio the 95% confidence radii. 





FlG. 9. — Histograms of the data (solid) and Model 4 (dashed). Model 4 predicts too few events at both small (p/a^ <, 0.25) and 
large (p/cr t ° ot ;> 3) separations. 
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Table 5 

Models Comparison Results for 4Br Events with Fluence Data 



Model 


Description 


Number of 
Parameters 


lo gio 
likelihood 


logio 
Odds Factor 





Minimal: Systematic 
Error Assumed to be 1.6° 





134.0 


134.0 


1 


Single Systematic Error 
applied to all locations 


1 


168.2 


166.9 


2 


Core-plus-Tail Systematic 
Errors 


3 


177.8 


175.0 


4 


Systematic Error a power 
law function of statistical 
uncertainty: 

^sy S =A<> stat /lT 


2 


170.4 


168.1 


8 


Systematic Error a power 
law function of fluence: 

0-sys=A((T stat /lxlO- 5 ) a 


2 


169.4 


166.8 



tively bright BATSE bursts (Fig. 2), it is difficult 
to use IPN data to test the locations of faint bursts. 
Graziani and Lamb (1996) used extrapolations of 
Models 4 and 8 to argue that faint BATSE bursts have 
large systematic errors, e.g., a sys = 8° for bursts with 
Cstat = 10°, and that the 3B catalog had larger sys- 
tematic errors for faint bursts than the IB catalog. 
Because our analysis finds no evidence that cr sys di- 
rectly depends on burst intensity, we consider such 
extrapolations unjustified. 

A more direct test of the location accuracy of weak 
bursts can be made using fluctuations from Cygnus 
X-l which triggered BATSE. These fluctuations are 
almost always near the trigger threshold. The spectra 
of these fluctuations are similar to GRBs in the en- 
ergy range (50-300 keV) used for determining loca- 
tions, but softer at higher energies. Because of the re- 
duced flux above the energy range used for burst lo- 
cations, the scattering model is less important for lo- 
cating Cygnus X-l than for GRBs; however, any in- 
tensity dependence of scattering other than the obvi- 
ous proportionality is implausible. Thirty-nine such 
events had an average statistical uncertainty <7 stat of 
12°, but an average separation from the true location 
of only 10°, implying a 95% confidence upper- limit 
ono- sys of7° (Meeganetal. 1996). 



There is, however, the possibility of a weak indi- 
rect correlation of cr sys with burst intensity. The frac- 
tion of events located using CONT data decreases 
with decreasing burst intensity (Fig. 10). The rea- 
sons for this trend are correlations between the opti- 
mum BATSE datatype for a location, GRB duration, 
and a stat : short events are best located with datatypes 
other than CONT to avoid adding background inter- 
vals to the burst data; additionally short events typ- 
ically have low fluences and thus larger values of 
(7 stat . In Model 33 the only difference between loca- 
tions obtained using CONT data and other locations 
is the value of the core systematic error. Coupled 
with the correlation between datatype and intensity, 
this datatype dependence of Model 33 implies that 
a larger fraction of weak bursts have a core system- 
atic error of 2.74° instead of 1.67°. This difference is 
unimportant for bursts with large values of cr stat . The 
systematic error is not a direct function of the statis- 
tical uncertainty and therefore does not extrapolate 
to large values for faint bursts. 

It seems plausible that the improvement of the 
models with direct intensity dependencies (Models 4 
and 8) over single-parameter models (0 or 1) is due 
to two correlations: the correlation between intensity 
and datatype used for the location and the probable 
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FlG. 10. — Fraction of events located using the CONT 
datatype versus <r stat . Each of the first 15 bins has 102 locations, 
while the last bin has 107. 
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correlation between datatype and systematic error. 
However, this may not be be the case since a model 
in which (j sys has no dependence on any burst prop- 
erty (Model 2) is favored over Model 4 by a large 
odds ratio. 

3.3. Previous BATSE Catalogs 

The locations of the previous BATSE catalogs, IB 
(Fishman et al. 1994), 2B (Meegan et al. 1994), 
3B (Meegan et al. 1995), and 4B (Paciesas et al. 
1997) are superseded by the revised locations of the 
4Br catalog (Paciesas et al. 1998), with the ex- 
ception that bursts of the 3B catalog which are not 
listed in the 4Br catalog paper have unaltered loca- 
tions which are considered to be 4Br locations. Key 
events affecting the quality of the locations were the 
deterioration of the data availability starting at the 
end of the IB catalog due to the failure of the tape 
recorders on-board CGRO, a substantial recovery of 
data availability at the end of the 2B catalog due to 
flight software improvements, and a major improve- 
ment in the location algorithm effective with the pub- 
lication of the 3B catalog (at which time all locations 
were revised). A lesser algorithmic improvement is 
implemented with the publication of the 4Br catalog 
(Paciesas et al. 1998), in which 208 locations are 
changed from the original 4B catalog. 

Table 6 presents analyses of these catalogs for his- 



torical interest and to aid in assessing work based 
upon these catalogs. Since the purpose is to under- 
stand BATSE location errors, in all cases the current 
IPN Supplement data are used rather than the pre- 
liminary IPN data available at the time of the pub- 
lication of a particular BATSE catalog. Only Mod- 
els 1 and 2 are shown because most of the catalogs 
are too small to support reliable analysis of models 
with more than about 3 parameters. The odds ra- 
tios favor Model 2 over Model 1 by factors ranging 
from 4 to 3 x 10 13 . The major algorithmic upgrade 
of LOCBURST can be seen in the comparison of the 
parameter values for the original 2B and 3B catalogs, 
or between the original 2B catalog and the equiva- 
lent subset of the 4Br catalog. The minor algorith- 
mic improvement between the 4B and 4Br catalogs 
is not apparent in Table 6, which is not surprising 
since only 13% of the locations were altered, mostly 
modestly. A more direct comparison validates the al- 
gorithmic upgrade used for the 4Br catalog (Paciesas 
et al. 1998). There is no evidence in the parameter 
values for subsets of the 4Br catalog for any time de- 
pendence of location quality, especially considering 
that the parameter values of Model 2 are correlated. 
In particular, the locations obtained during the pe- 
riod of lesser-quality data (2B-1B) have error model 
parameters consistent with the other time periods. 

4. CONCLUSIONS 
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Table 6 

Error Models for Various Datasets 

Num. of Num. Num. Model 1 Model 2 

Trigger BATSE w. IPN with a sys fi a* ys 

Catalog Range GRBs data 7 (deg) (deg) (deg) 



Superseded catalogs analyzed using originally published data: 



IB 105 1466 260 54 6 3.4 ± 0.4 1.6 ±0.4 0.57^ ^ 5.0^ 

2B 105 2230 585 128 18 4.5 ± 0.3 1.8 ±0.3 0.75 ± 0.07 9 3+H 

3B 105 3174 1122 264 19 2.88 ±0.15 1.69 ±0.22 0.71 ±0.11 4.9 ± 0.8 

4B 105 5586 1637 411 19 2.90 ±0.12 1.84 ±0.18 0.74 ±0.09 5.0 ± 0.7 



Subsets of the 4Br catalog: 



full 


105 


5586 


1637 


411 


19 


2.80 ±0.12 


1.85 ±0.16 


0.78 ±0.08 


r 1+0.8 
J- 1 -0.6 


IB 


105 


1466 


263 


54 


6 


2.56 ± 0.30 


1.6±0.3 


7Q+ 010 


c 4+2.5 


2B 


105 


2230 


586 


128 


18 


2.80 ±0.20 


1.55 ±0.26 


n 79+0.09 

u./z_ 13 


5 2 +L2 


2B-1B 


1467 


2230 


323 


74 


12 


3.0 ±0.3 


1.5 ±0.4 


0.67 ±0.18 


5 1 +L5 
J - 1 -i.o 


3B 


105 


3174 


1122 


264 


19 


2.89 ±0.15 


1.77 ±0.18 


0.76 ± 0.09 


c 4+O.9 


4B-2B 


2232 


5586 


1051 


283 


1 


2.80 ±0.14 


2.08 ± 0.21 


0.843E 


5 4 +L3 


4B-3B 


3177 


5586 


516 


147 





2.61 ± 0.20 


2.1 ±0.4 


0.81!^ 


4 4+3.8 



This paper presents an improved model of the 
BATSE GRB location errors. Future work may 
incorporate improvements such as the incorpora- 
tion of more locations (BATSE and IPN), addi- 
tional 'point' IPN locations when more interplan- 
etary spacecraft, such as the Near Earth Asteroid 
Rendezvous (NEAR) mission, are added to the net- 
work, the use of locations obtained with BeppoSAX 
(which was launched several months before the end 
of the 4Br catalog) and HETE II, and the extension of 
the analysis to use comparison locations with errors 
not negligible compared to BATSE (e.g., COMPTEL 
and WATCH locations). Another possible extension 
would be to use elliptical error boxes rather than the 
azimuthally symmetric distributions assumed here. 
We continue to search for correlations between lo- 
cation errors and other properties. Such correlations 
might explain why some bursts are in the tail of the 
location distribution. If we identify a correlation, we 
might be able to use that knowledge to improve the 
location algorithm. 

The parameter values of model fits to catalogs pre- 
ceding the 4Br catalog confirm that the location al- 
gorithm introduced with the 3B catalog was a signif- 
icant improvement over the previous algorithm. 

The models of this paper are highly significant 



improvements over previous models of the BATSE 
GRB location error distribution. We recommend the 
three-parameter 'core-plus-tail' model (Model 2) for 
most analyses. There is evidence that a 'core-plus- 
tail' model with core systematic error depending on 
datatype (Model 33) is superior, and this model can 
be used for the most critical applications. 

A promising use of BATSE and IPN locations is 
rapid (minutes to hours) optical observations of GRB 
locations obtained with the BATSE Rapid Burst Re- 
sponse (RBR) System (Kippen et al. 1998), fol- 
lowed by detailed analysis of the smaller error box 
obtained either with the intersection of a single IPN 
annulus with the BATSE box or the very small er- 
ror box obtained by the intersection of a pair of IPN 
annuli. The RBR locations are obtained by human 
operators using the standard LOCBURST algorithm 
and have accuracies consistent with the accuracies of 
BATSE catalog locations. 

MSB, GNP, RMK, and JJB acknowledge NASA 
Grant NCC8-65 for support. KH is grateful to JPL 
Contract 958056 for support of Ulysses operations 
and to NASA Grant NAG5-1560 for support of the 
IPN. We appreciate the suggestions of the anony- 
mous referee. 
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APPENDIX 



DERIVATION OF P(p) 

The goal is to test the probability density function p(^) for the separation 7 between the BATSE location and 
the true location. In most cases only single IPN annuli are available and therefore only the smallest distances 
p between the BATSE location and the IPN annuli can be determined. In such cases the testing of p(j) must 
be done indirectly by comparing p(p) to the observations. The subject of this Appendix is deriving p(p) given 
p("f), assumed to be the Fisher distribution (eq. 5). 

Figure 11 depicts the geometry of comparing a BATSE location to a single IPN annulus. All circles and 
positions are on the surface of the unit sphere. The IPN annulus, which may not be a great circle, has radius 
r and is assumed to have a negligible width Sr. The coordinate system is chosen so that the center of the IPN 
annulus is at z and the true location t is in the xz plane. The BATSE location B is at an arbitrary position on 
the unit sphere. 

Locations may be expressed in spherical coordinates (9 7 <f>), where 9 is the angle from z and <f> is the azimuthal 
angle from x (see Fig. 1 1). Since 

p=\9-r\, (Al) 

the probability density function p(p) can be readily obtained from p(9). The assumed probability density 
function is a function of 7, so it is necessary to derive the the relationship between the the coordinate system 
(6>, 4>) and the system (7, ip), where 7 is the angle from T to B and ip is an azimuthal angle about T. 

f = sinri+cosrz, (A2) 



B = sin # cos i+ sin 6> sin + cos 6>z. (A3) 

Therefore, 

COS7 = TB (A4) 

= cos(6'-r)-sin^sinr(l -cos0). (A5) 

Since the transformation is a rotation, the Jacobean determinant is unity and the transformation between 
probability density functions is simply p{0,4>) = ^(7,^). The probability density function is assumed to be 
rotationally symmetric, so pdy,ip) = p(7). 

Substituting for cos 7 from eq. A5 into eq. 5 for the Fisher distribution, 

p(9, (f)) dQ = e «co S (0-r) e - K sin0sinr(l-cos0) JQ (A6) 

2n(e K -e K ) 




Fig. All. — Diagram of the geometry relating the true lo- 
cation T to the BATSE location B. The coordinate system is 
chosen so that the center of the IPN annulus is at z and the true 
location t is in the xz -plane. The annulus has radius r. The 
BATSE location B is at spherical coordinates 9, <\>. The angle 
between the true location t and the BATSE location B is 7, 
while the angle between the BATSE location B and the closest 
location on the annulus C is p. 
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and 



p(9) d9 = f V d(f> p(9, (j)) d6 (A7) 
Jo 

/"rf0e- Ksinesinr(1 - cos * ) dd. (A8) 
Jo 



_^ KCOs(6»-r) 



7r(e K -e K ) 

Finally, using formula 8.431.3 of Gradshteyn and Ryzhik (1980), 

m de = e-^) /Q(KSin f Sinr) dO. (A9) 

r qK g— k qK sin a sin r 

Including the factor Q Ksm0smr in the last term makes that term more amenable to numerical evaluation (e.g., 
routine BESSIO of Press et al. (1992) is easily modified in this manner). 

Usually there are two values of 9 for a value of p, corresponding to locations inside and outside of the 
annulus, so that 

p(p) dp = p(9 = r-p) d9+p(9 = r+p)d9. (A 10) 

In the unusual case that p > r, only the 9 = r + p term contributes. 
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